Rare K Decays: Results and Prospects 



Laurence Littenberg 

Brookhaven National Laboratory, Upton, NY 11973 

Abstract. Recent results on rare kaon decays are reviewed and prospects for future experiments are 
discussed. 



INTRODUCTION 

In recent years the study of the rare decays of kaons has had three primary motivations. 
The first is the search for physics beyond the Standard Model (BSM). Virtually all at- 
tempts to redress the theoretical shortcomings of the Standard Model (SM) predict some 
degree of lepton flavor violation (LFV). Decays such as Kl /j^e^ have very good 
experimental signatures and can consequently be pursued to remarkable sensitivities. 
These sensitivities correspond to extremely high energy scales in models where the only 
suppression is that of the mass of the exchanged field. There are also theories that pre- 
dict new particles created in kaon decay or the violation of symmetries other than lepton 
flavor. 

The second is the potential of decays that are allowed but that are extremely sup- 
pressed in the SM. In several of these, the leading component is a G.I.M.-suppressed[[l|] 
one-loop process that is quite sensitive to fundamental SM parameters such as Vtd- These 
decays are also potentially very sensitive to BSM physics. 

Finally there are a number of long-distance-dominated decays which can test the- 
oretical techniques such as chiral Lagrangians that purport to explain the low-energy 
behavior of QCD. Knowledge of some of these decays is also needed to extract more 
fundamental information from certain of the one-loop processes. 

This field is quite active as indicated by Table |T] that lists the decays for which results 
have been forthcoming in the last couple of years as well as those that are under analysis. 
Thus in a short review such as this, one must be quite selective. 



BEYOND THE STANDARD MODEL 

There were several K decay experiments dedicated to lepton flavor violation at the 
Brookhaven AGS during the 1990's. These advanced the sensitivity to such processes by 
many orders of magnitude. In addition, several "by-product" results on LFV and other 
BSM topics have emerged from the other kaon decay experiments of this period. Table ^ 
summarizes the status of work on BSM probes in kaon decay. 



TABLE 1. Rare K decay modes under recent or on-going study. 
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TABLE 2. Current 90% CL limits on K decay modes violating the SM. The 
violation codes are "LF" for lepton flavor, "LN" for lepton number, "G" for 
generation number, "H" for helicity, "N" requires new particle 
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It is clear from this table that any deviation from the SM must be highly suppressed. 
The LFV probes in particular have become the victims of their own success. The specific 
theories they were designed to test have been killed or at least forced to retreat to the 
point where meaningful tests in the kaon system would be very difficult. Both kaon 
flux and rejection of background are becoming problematical. Analysis of data already 
collected is continuing but no new kaon experiments focus sed on LFV are being planned. 
Interest in probing LFV has migrated to the muon sector. 



ONE LOOP DECAYS 

In the kaon sector experimental effort has shifted from LFV to "one-loop" decays. 
These are GIM-suppressed decays in which loops containing weak bosons and heavy 
quarks dominate or at least contribute measurably to the rate. These processes include 
Kl Ti^vv, TT+vv, Kl — > fj^fj^, Kl — > n^e^e^ and Kl ttV^^u". In some cases 

the one-loop contributions violate CP. In one, Kl n^vv, this contribution completely 
dominates the decay [jl^. Since the GIM-mechanism tends to enhance the contribution 
of heavy quarks in the loops, in the SM these decays are sensitive to the product of cou- 
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FIGURE 1. K decays and the unitarity plane. The usual unitarity triangle is dashed. The triangle that 
can be constructed from rare K decays is solid. See text for further details. 



plings Vi*Vtd, often abbreviated as Xt. Although one can readily analyze these decays in 
terms of the real and imaginary parts of Xt, for comparison with results in the B system, 
it is conventional to parameterize them in terms of the Wolfenstein variables, A, p, and 
r\. Fig. [H shows the relationship of rare kaon decays to the unitarity triangle construc- 
tion. The dashed triangle is the usual one derived from V*^Vud + ^cb^cd + ^tb^td = 0, 
whereas the solid triangle illustrates the information available from rare kaon decays. 
Note that the "unitarity point" at the apex, (p,Ti), can be determined from either tri- 
angle, and disagreement between the K and B determinations implies physics beyond 
the SM. In Fig. |T]the branching ratio closest to each side of the triangle determines the 
length of that side. The arrows leading outward from those branching ratios point to 
processes that need to be studied either because they potentially constitute backgrounds, 
or because knowledge of them is required to relate the innermost branching ratios to 
fundamental parameters. Ki /J^/J , which can determine the bottom of the triangle 
(p), is the process for which the experimental data is the best, but for which the the- 
ory is most problematical. — > tt^vv, which determines the height of the triangle is 
theoretically the cleanest, but experiment is many orders of magnitude short of the SM- 
predicted level. 7r+vv, which determines the hypotenuse, is nearly as clean as 

Ki TT^vv and has been observed. Prospects for tt+vv are probably the best of 

the three since it is already clear it can be exploited. 



The short distance component of this decay can be quite reliably calculated in the 
SM[|TT|]. The most recent measurement of its branching ratio [[T2|] based on some 6200 



events gave B{Ki /j^/j ) = (7. 18 ±0. 17) x 10 ^. However is dominated 

by long distance effects, the largest of which, the absorptive contribution mediated by 
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FIGURE 2. Components of the calculation of . 



Kl 77, accounts for (7.07 ±0. 18) x 10"^. Subtracting the two, yields a 90% CL upper 
limit on the total dispersive part of B(Ki /J^/J ) of 0.37 x 10^. One can do a little 
better than this in the following way. The actual quantity measured in Ref [ [T^ ] was 

B{Ki-<-n+Tc-) ~ (3-48 ±0.05) X 10^^ One wants to subtract from this measured quantity 

the ratio . Fig ^ shows the components of this latter ratio, obtained from 

Ref. [0, whose product is (3.435 ±0.065) x lO"*^. 

The subtraction yields '^'s'^l'^+I^')^ = (0-045 ±0.082) x 10"^ (where B'^'^p refers to 
the dispersive part of B{Kl — > ir^^)). g(^^^Z^T^+Tf^ can then be multiplied by B{Ki 

71+71-) = (2.056 ±0.033) X lO^^to obtain 5^"'' (iSC^^A'+A'") = (0.093 ±0.169) x lO^^, 
or B'^'^P{Ki fJ^fJ') < 0.31 X 10"^ at 90% CL. Note that some of the components 
represent quite old measurements. Since B{Ki — > ^^^^) and B^^{Kl fJ^fJ^) are so 
close, small shifts in the component values could have relatively large consequences 
for B^'^P{Ki — i- Several of the components could be remeasured by experiments 

presently in progress]]. Now if one inserts the result of even very conservative recent 
CKM fits into the formula for the short distance part of B{Ki fj~^fj~), one gets 
rather poor agreement with the limit of B'^''^p{Ki /J^/J^) derived above. For example 
the 95% CL fit of Hocker et al.[|l5l][|l6l], p = 0.07-0.37, gives B^^{Kt ^i+n') = 
(0.4 — 1 .3) X 10^^. So why haven't we been hearing about this apparent violation of the 
SM? 

The answer is that unfortunately Kl 7*7* also gives rise to a dispersive contribution, 
which is much less tractable than the absorptive part, and which can interfere with the 
short-distance weak contribution that one is trying to extract. The problem in calculating 
this contribution is the necessity of including intermediate states with virtual photons of 



There is a new preliminary result from the KLOE experiment of B{Ks %^%^)/B{Ks K^K^) = 
2.192 ± 0.003stat ± 0.016^v,st []l4|]. Fortuitously, inserting this result in place of the PDG value makes no 
difference to the final limit. 
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FIGURE 3. Left: Spectrum of x = {m^^/mKY in Ki /j'^/j yfrom Ref.[p7|. Right: Determinations of 
the BMS parameter Uk* from three Kl decays involving virtual photons. 



all effective masses. Thus such calculations can only partially be validated by studies 
of processes containing virtual photons in the final state. Recently there have been 
publications on 



YA^V^O (9327 events), Kl e+e^e^e^]^] (441 events), and 
(38 events) and there exist slightly older high statistics data on 



Kl ye^e [ [T9| ] (6864 events). Figure |^-left shows the spectrum of x = {mjj^/mK) 



from Ref.QTTl]. The disagreement between the data (filled circles with error bars) and 
the prediction of pointlike behavior (histogram) clearly indicates the presence of a form 
factor. A long-standing candidate for this is provided by the BMS model [|20|] which 
depends on a single parameter, a^:* . 

Fig. fright shows three determinations of this parameter. The level of agreement of 
these results leaves something to be desired. Fitting to a more recent parameterization 
of these decays[|2Tl] also results in quite marginal agreement. This may improve when 
radiative corrections are properly taken into account. Thus additional effort, both exper- 
imental and theoretical, is required before the quite precise data on B{Kl l-i^l-i ) can 
be fully exploited. 
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Theoretically K^ tt+vv is remarkably clean, suffering from none of the long dis- 
tance complications of Kl — »• /J~^/J ■ The hadronic matrix element, so often a prob- 
lem in other processes, can be calculated to a ~ 2% via an isospin transformation 



from that of ^e3[E3]. Interest in K^ tt+vv is driven in large part by its sensi- 



tivity to Vfd (it is actually directly sensitive to the quantity |yf*y;rf|). Its amplitude is 
proportional to the dark slanted line at the right in Fig. |l]. This is equal to the vec- 
tor sum of the line proportional to iVf^/l/A^i^ (where X = sinQcabibho) and that from 
(1,0) to the point marked po. The length po — 1 along the real axis is proportional to 
the amplitude for the charm contribution to K^ tt+vv. The QCD corrections to this 
amplitude, which are responsible for the largest uncertainty in B{K^ 71+ vv), have 
been calculated to NLLA[jTT|]. The residual uncertainty in the charm amplitude is es- 
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FIGURE 4. Left: new K'^ 7t+vv event. Right: Range vs energy of 7t+ in the final sample. The circles 
are 1998 data and the triangles 1995-7 data. The events around £ = 108 MeV are — > background. 
The simulated distribution of expected signal events is indicated by dots. 



timated to be ~ 15% which leads to only a ~ 6% uncertainty p3|] in extracting \Vtd\ 
from B{K^ tt+vv). Recently AGS E787 has seen evidence for a second event of 
TT+vv [0] (see Fig. which, combined with previous data p^], yields a branch- 
ing ratio B{K^ 7r+vv)= (1.57+Q g2) x 10"^*^. By comparison, a fit to the CKM phe- 
nomenology yields the expectation (0.72 ±0.21) x 10"^"[^. It is notable that E787 
has established methods to reduce the residual background to ~ 10% of the the signal 
branching ratio predicted by the SM. 

A new experiment, AGS E949[^, based on an upgrade of the E787 detector, is about 
to begin its first physics run. Using the entire flux of the AGS for 6000 hours, it is 
designed to reach a sensitivity of ~ lO^^Vevent. In June 2001, Fermilab gave Stage 1 
approval to an experiment (CKM[^7|]) to extend the study of TT+vv by yet an 

another order of magnitude in sensitivity. This experiment, unlike all previous ones on 
this process, uses an in-flight rather than a stopping technique. This experiment is 
expected to start collecting data in 2007 or 2008. 

Fig. ^ shows the history and expectations of progress in studying tt+vv. 



Ki 71"VV 

Ki — > tt'^VV is the most attractive target in the kaon system, since it is direct CP- 
violating to a very good approximation [[TOj |28| ] (B{Ki — > 7r*^vv) <^ T]^). Like 
TT+vv it has a hadronic matrix element that can be obtained from Kgi,, but, it has no 
significant contribution from charm. As a result, the intrinsic theoretical uncertainty 
connecting B{Ki tt^vv) to the fundamental SM parameters is only about 2%. Note 
also that B{Ki tt^vv) is directly proportional to the square of ImXt and that ImXt = 
—J/ — -y)] where J7 is the Jarlskog invariant[p9|]. Thus a measurement of B{Kl — ^ 
7r*^vv) determines the area of the unitarity triangles with a precision twice as good as 
that on B{Kl 7l°vv) itself. 
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FIGURE 5. History and prospects for the study of 7t+vv. Points without error bars are single 

event sensitivities, those with error bars are measured branching ratio. 



B{K]^ — > 7r*^vv) can be bounded indirectly by measurements of B{K^ 7r''"vv) 
through a nearly model-independent relationship pointed out by Grossman and Nir[|30[]. 
The application of this to the new E787 resuh yields B{Ki^ 7r°vv)< 1.7 x lO^at 
90% CL. This is far tighter than the current direct experimental limit, 5.9 x 10^, ob- 
tained by KTeV[0. To actually measure B{Kl TT^Vv) at the SM level (~ 3 x 10'^^), 
one will need to improve on this by some five orders of magnitude. The KEK E391a 
experiment proposes to achieve a sensitivity of ~ 3 x 10^'*^/event which would bet- 
ter the indirect limit by a factor five, but would not quite bridge this gap. It will serve 
as a test for a future much more sensitive experiment to be performed at the Japanese 
Hadron Facility. E391a features a carefully designed "pencil" beam, and a very high 
performance photon veto. The active photon detector is a CsTpure crystal calorimeter. 
The entire rather compact apparatus will operate in vacuum. Beamline construction and 
tuning started in March 2000 and physics running is expected to begin in Fall, 2003. 

The KOPIO experiment at BNL (E926) takes a completely different approach, 
exploiting the intensity and flexibility of the AGS to make a high-flux, low-energy, 
microbunched Ki beam. The proposed experiment is shown in Fig. ^. The neutral beam 
will be extracted at ~ 45° to soften the spectrum sufficiently to permit time-of- 
flight determination of the velocity. The large production angle also softens the 
neutron spectrum so that they (and the K^) are by and large below threshold for the 
hadro-production of TT^'s. The beam region will be evacuated to 10^ Torr to further 
minimize such production. With a 10m beam channel and this low energy beam, the 
contribution of hyperons to the background will be negligible. Ki decays from a ~ 3m 
fiducial region will be accepted. Signal photons impinge on a 2 Xq thick preradiator 



capable of measuring their direction. An alternating drift chamber/scintillator plane 
structure will allow energy measurement as well. A high-precision shashlyk calorimeter 
downstream of the preradiator will complete the energy measurement. The photon 
directional information will allow the decay vertex position to be determined. Combined 
with the target position and time of flight information, this provides a measurement of 
the Kt 3 -momentum so that kinematic constraints as well as photon vetoing are available 
to suppress backgrounds. The leading expected background is Kt n^n^, which is 
initially some eight orders of magnitude larger than the predicted signal. However since 
7r*^'s from this background have a unique energy in the Ki center of mass, a very 
effective kinematic cut can be applied. This reduces the load on the photon veto system 
surrounding the decay region to the point where the hermetic veto techniques proven 
in E787 are sufficient. In fact most of the techniques necessary for KOPIO have been 
proven in previous experiments or in prototype tests. KOPIO aims to collect about 50 
Ki — > 7r*^vv events with a signal to background ratio of 2:1. This will permit T] to be 
determined to ~ 10%, given expected progress in measuring nit and Vcb- KOPIO will 
run during the ~20 hours/day the AGS is not needed for injection into RHIC. 
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FIGURE 6. Layout of the KOPIO detector. 



These are reactions initially thought experimentally more tractable than Ki n^\\. 
Like Kl n^v, in the SM they are sensitive to ImXt, but in general they have 
different sensitivity to BSM effects Q]. Although their signatures are intrinsi- 
cally superior to that of — > 7l*^vv, they are subject to a serious background 
that has no analogue in the case of the latter: Ki — » 7y£+£^. This process, a ra- 
diative correction to Ki yi^i^, occurs roughly 10^ times more frequently than 
Kl n^i^i^ . Kinematic cuts are quite effective, but it is very difficult to improve 



the signakbackground beyond about 1 : 5[|35|]. Both varieties of Ki flVl have 



been observed, B{Kl W+(^^)ky>5MeV = {5M±0.15{stat) ±0.32{syst)) x 10"^[|3|] 
and B{Kl ^ m^l^-),nr,>iMeV/c2 = {lOAt]i{stat)±OJ{syst)) x IQ-'^^; both 
agree with theoretical prediction. By comparison, in the SM B'^"''^'^' (K^ TT^e+e^) is 
predicted to be ^] (4.3 ± 2.1) x 10" ^nd B'^'''''\Kl n^p+p-) about fives times 
smaller. 

In addition to this background, there are two other issues that make the extraction 
of short-distance information from Ki n'^i'^i^ problematical. First, there is an 
indirect CP-violating contribution from the component of K^ given by \e\^'^B(Ks 

vPe^e^) which is of the same order of magnitude as the direct CP- violating piece^. 
The exact size of this contribution will be predictable if and when B{Ks n^e^e^) is 



measured, hopefully by the upcoming NA48/1 experiment [|39|]. Second is yet another 
contribution of similar size mediated by Kl — > ;r^yy which is CP-conserving. To some 
extent this contribution can be predicted from measurements of the branching ratio 
and kinematic distributions of Ki 7i°7y, and thousands of these events have been 
observed. However as indicated in Table ^ a new result from NA48[pO|] disagrees by 



nearly 3a from the previous result from KTeVfPT]]. The change in the vector meson 
exchange contribution, characterized by the parameter ay, reduces the predicted size 
of B*^^^"'"^ {Kl Ti^e^e^) considerably which is good news for the prospects of 
measuring B'^"''^^'' {Ki — > n'^e^e^). However the validity of the current technique for 
predicting B'"^^'^""^{Ki — > n^e^e^) from Ki 7r*^yy has recently been reexamined [|43|], 
and questions raised about the functions used to fit the spectrum and about the treatment 
of the dispersive contribution. Thus both the theoretical and experimental situations are 
quite unsettled at the moment. Depending on whose data and whose theory one uses, 
values from 0.25 x 10"'^ 7 3 ^ 10-12 are predicted for B^^-''""' {Kl n^e+e-). 



TABLE 3. Results on Kl Ti^yy. 



Exp. B{Kl TTTf) X 10^ av Ref. 



KTeV 1.68±0.07,,„,±0.08,y,, -0.72 ±0.05 ±0.06 

NA48 1.36±0.03,,„,±0.03,y,,±0.03„„„„ -0.46 ± 0.03 ± 0.03 ±0.02,w |4C] 



The current experimental status of Kl n^E^i is summarized in Table A factor 
2.5 more data is expected from the KTeV 1999 run, but as can be seen from the table. 



There is also an interference term between the indirect and direct CP-violating amplitudes. 
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FIGURE 7. /j+ polarizations in % /j^/j , plotted against the muon cm energies. Left: longitudinal 
polarization. Right: out-of-plane polarization. 



background is already starting to be observed at a sensitivity roughly 100 times short of 
the expected signal level. 

TABLE 4. Results on Kl K°£+£-. 
Mode 90% CL upper limit Est. bkgnd. Obs. evts. Ref. 



Ki^n^e+e- 5.1x10-1° 1.06 ±0.41 2 Q 

Ki^n^fj+fj- 3.8 xlO-"* 0.87 ±0.15 2 g| 

To make a useful measurement under these conditions will require markedly increased 
statistics on the signal and both theoretical and experimental advances in the ancillary 
modes n^YY Ks — > TT^e+e^. Various approaches for mitigating these prob- 

lems have been suggested over the years including studies of the Dalitz Plot [Q, the 
time development p7|], or both [Q. However an innovative approach has recently been 
suggested in which muon polarization in Ki — > ttV^A' well as kinematic dis- 
tributions are exploited. The /u^ longitudinal polarization is proportional to the direct 
CP-violating amplitude, whereas the energy asymmetry and the out-of-plane polariza- 
tion depend on both indirect and direct CP violating amplitudes. As shown in Fig. 0, 
the polarizations involved turn out to be extremely large so that enormous numbers of 
events may not be required. 



An alternative parameterization 

Although it is customary to write the branching ratios and other observables of the 
one-loop processes in terms of the Wolfenstein parameterization of the CKM matrix, 
this parameterization is not really natural to the kaon system, and puts results from this 
system at a certain disadvantage in comparisons with those from the B system. To extract 
information on p and r\, for example, it is necessary to divide the physical measurements 
by X^A^, thereby introducing "external" contributions to the uncertainty of Sox, and 4oa- 
One can avoid this by resorting to expressions for the branching ratios in terms of the 



quantity Xf Since as noted above, the imaginary part of this quantity determines the area 
of all unitarity triangles, it is no less fundamental than p and Tj. 

The formulae for the branching ratios of three of the decays discussed above are: 



where 



5(A-+^7C+vv) = ^{XcX+Re(kt)Xtf + {.Im{Xt)Xtf] (1) 
B{Ki^^Tl%v) = ^'{Im{h)X,f (3) 



w 



1.55 X lO"'^ 
6.32 X 10^3 



X/^+y„227r25m40 



6.77 X 10 



-4 



(4) 
(5) 
(6) 



and to a good approximation the Inami-Lim [ |50| ] functions Xt = 1.56(mf/170GeV)^-^^, 
Yt = 1.02(mf/170GeV)i-^^. The quantities = 0.901 and r/^^ = 0.944 are isospin 
correction factors that relate the hadronic matrix elements of the K ttvv processes 
to that of — > n^e^v . The terms X and Y/v/^ are the Inami-Lim functions for the 



charm contributions which, after correction to NLLA, are known to about 15%. 
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FIGURE 8. Comparison of 90% CL constraints from current data on rare kaon decays with 95% CL 
constraints from a typical unitarity fit (based on Ref. JT^). The allowed region from kaon decays lies 
between the two circles and within the outer two vertical lines. 

Fig. H shows the 90% CL constraints currently available from l^^l^^ and 

To extract a limit from the former we adopt the value for the maximum 
long distance dispersive contribution from Ref. [|2l]]. Also shown is the region in the 



X/ plane bounded by the CKM fit mentioned above [|T^. The two kinds of information 
are clearly consistent at the moment. Note, however, that if the current central value 



of B{K^ TT+VV) should hold through E949, and expected progress is made in the B 
sector, this agreement could prove short-lived, as shown on the left of Fig. When, 
eventually, we have 10-15% measurements of B{K^ tt+vv) and B{Kl — > Ti^vv), 
comparison of K and B results will become a critical test of the SM. Fig. ^ (right) 
illustrates a scenario in which such a failure is evident. 
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FIGURE 9. Left: Similar plot to Fig. || after 10^'Vevent 7t+vv experiment. Assumes central 

value of B{K^ 7i+vv) stays the same and also that CKM fit contours and precision on m, are improved 
by a factor 2. Right: Similar plot for possible scenario after 10% measurements of | A,, | and Im{Xt). Further 
improvements in CKM parameters and assumed. 



CONCLUSIONS 

The success of lepton flavor violation experiments in reaching sensitivities correspond- 
ing to mass scales of well over 100 TeV has helped kill most models predicting acces- 
sible LFV in kaon decay. The most popular varieties of SUSY predict LFV at levels far 
beyond the current experimental state of the art []5T|]. Thus new dedicated experiments 
in this area are unlikely in the near future. 

The existing precision measurement of /J^/J will be very useful if theorists can 
make enough progress on calculating the dispersive long-distance amplitude, perhaps 
helped by experimental progress in yi^i^, A leptons, etc. The exploitation 
of Ki — > /j^/j would also be aided by higher precision measurements of some of the 
normalizing reactions, such as yy. 

— > TT+vv will clearly be further exploited. Two coordinated initiatives are de- 
voted to this: a lO^^Vevent experiment (E949) just underway at the BNL AGS and a 
10^^^/event experiment (CKM) recently approved for the FNAL Main Injector. The first 
dedicated experiment to seek Ki tt^vv (E391a) is proceeding and an experiment (KO- 
PIO) at the AGS with the goal of making a ~ 10% measurement of Im(kt) is approved 
and in R&D. 

Measurements of tt+VV and Ki tt^VV can determine an alternative unitarity 

triangle that will offer a critical comparison with results from the B system. If new 
physics is in play in the flavor sector, the two triangles will almost certainly disagree. 
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